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Abstract. In this paper a new optimum tuning method of PI controllers in 
first-order time-delay systems, based on the deadbeat response to a step set- 
point variation, is presented. The deadbeat performance, already studied for 
the plants without delay, consists of a fast achievement of the steady state with 
an overshoot included in a narrow band and with minimum rise and settling 
times. In the proposed method the rise and settling times are both replaced 
with the integral of the squared error and the constraint, which requires a 
controller output less than a preset value, is added. The transient behaviour 
is strictly evaluated by means of the analytical solutions of the relative differ- 
ential difference equations, obtained by the method of steps for the first time. 
Moreover the proposed tuning method is compared with other three methods, 
selected among the most used, which are the Ziegler-Nichols time-domain and 
frequency-domain and the Zhuang and Atherton ISTE ones. Finally tuning 
charts, having as coordinates the two PI controller parameters and provided 
with the borderlines of the stability region, are introduced. 



1. Introduction 

The implementation of an optimum setting of the proportional integral derivative 
(PID) controllers , especially in time-delay systems, is one of the most attractive 
goal of the control designers and it has been investigated from more than 60 years. It 
seems unbelievable that today an original and interesting word can be said, but this 
happens in the present paper for the PI controllers thanks to a new mathematical 
tool, which is the explicit solution of differential difference equations detailed in pQ . 

Before the enunciation of the proposed method, let us mention some of the most 
significative existing ones. The first attempt has been performed by Ziegler and 
Nichols, who presented two design methods, one time-domain and one frequency- 
domain 2 J. In both methods the amplitude of the decay ratio, which is the ratio 
between two consecutive maxima of the control error during a load disturbance, 
has been fixed equal to one quarter. The Ziegler-Nichols tuning rules have the ad- 
vantage of being easy to use, but give closed-loop systems with very poor damping. 
Modified expressions of the controller parameters, which improve the damping, 
have been proposed by Chien, Hrones and Reswick [3] and by Cohen and Coon 
A second method, named Kappa- Tau [5], assumes the maximum sensivity as 
design parameter and applies the dominant pole design to typical process control 
models. A third method is based on the minimization of an integral performance 
criterion, whose the most well known is the integral of the squared error (ISE). 
Zhuan and Atherton [6] have approximated the exponential term of the first-order 
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time-delay plant model using a 3/3 order Pade and have minimized the ISE, ISTE 
and IST 2 E performance criteria for setpoint and load changes. An other method, 
which deserves to be cited, is the deadbeat response one, studied for processes 
without time delay and defined in para. fO.2.5 of [7] pag. 171 as follows: 

• Zero steady-state error 

• Minimum rise time and settling time 

• 0.1% <= controlled variable overshoot <= 2% 

• Controlled variable undershoot <= 2% 

The proposed method, which is a modification of the deadbeat response one, can 
be defined concisely, but exhaustively, as follows: it minimizes the integral squared 
error (ISE), during a closed- loop step setpoint change, with the controlled variable 
overshoot (PO y ) equal to a preset value (for example 1.05%) and the controller 
output overshoot (PO v ) lower than a preset value (for example 10%). Only two 
performances indexes, ISE and PO y or PO v , are involved in this tuning optimiza- 
tion and therefore a PI controller, provided with two parameters, is adequate i.e. a 
PID controller is not required. 

A significative advantage of this method is the avoidance of the controller windup 
without any additional sophisticated scheme. It is known that one of the nonlinear 
issues of any controller is the actuator saturation, since from the practical point 
of view the actuator transfer function is bonded and the controller output may be 
out of the rated band. When the actuator saturates, the feedback loop is broken 
and the actuator output is constant and then independent of the error signal. In 
this case a large overshoot of the controlled variable may occur, a long time may be 
needed for the recovery of the system and the actuator may bounce several times 
between the two bond values. If the rated value of the controller output is set lower 
than the extreme point of the linear part of the actuator transfer function and if 
the controller is tuned in such a way that the known maximum overshoot of the 
controller output never reaches this extreme point, none windup occurs. 

The feedback system taken into consideration is represented by the block diagram 
of Fig. 1 and the transfer functions are mathematically described by 

(i.i) G ^= K ^L 



(1.2) G c (s) = K p + 



Ki 
s 



(1.3) G c (s) = K p (l + i- 

where K is the plant steady-state gain, T p the positive plant time constant, L the 
plant time delay and K p , Ki and Ti the parameters of the PI controller. 

In order to get equations independent of the real values of the parameters, the 
normalized time t and dimcnsionlcss parameters are introduced as follows: 

• Systems without delay 

t referred to the plant time costant T p and h = K K p , ti = Ti/T p . 

• Systems with delay 

t referred to the plant time delay L and t p = T p / L, h = K K p , hi — K KiL. 
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Figure 1 . Feedback control system 



The open loop transfer functions for this system without time delay F% (a) (a = 
T p s in (fTTTj) and (|L3]0 and with time delay F 2 (a) {a = L s in (fl~l~j) and (fO|) ) are 
given respectively by 



(1.4) 



Jk(tr) 



/i(l +ti<j) 
Ua{l + a) 



(1.5) 



h a 



a (1 + ipcr) 



The tuning points will be positioned in tuning charts, whose coordinates are h 
and ti for systems without delay and h and hi for systems with delay. The stability 
region borderlines and the curves of constant phase margin will be also plotted in 
these charts and therefore the stability of each selected tuning point may be easily 
verified. 

In the next Sections the proposed tuning method will be detailed, applied to 
both systems, without and with time delay, and finally compared with the two 
Ziegler-Nichols and with the Zhuang-Atherton ISTE for setpoint change ones. 



2. Systems without time delay 

For h+ 1 > the system is stable and hence all the tuning points of the quadrant 
h > and ti > are stable. The equations of the curves T y , T v and T i: relative 
respectively to the controlled variable overshoot PO y , the controller output over- 
shoot PO v and the integral of the square error ISE, are given by (|A.7|) , (|A.8|) and 
(|A.9jl . Equating (|A.3|) to zero, one obtains for the damping borderline 

Ah 

U = 



(i + hy 



The response is overdamped if the selected point is above this line and underdamped 
if below. In Fig. 2 the following curves are plotted: 

• r^: damping borderline 

• I\,: PO v = 0.1 

• T y2 , T y3 : respectively PO y = 0.1%, PO y = 1.05%, PO y = 2% 

• T a , r. i2 , r i3 : respectively ISE = 1.0, ISE = 1.2, ISE = 1.4. 

The optimum tuning point, named B and corresponding to the minimum value of 
the ISE for PO y — 1.05% and PO v <= 0.1, is the intersection between T y 2 and 
T,;, since the ISE always decreases when PO v increases for a given value of PO y . 
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T ( i : damping bor der line; T v : P0 V = 0.10 

T yl : PO y = 0.001; T y2 : PO y = 0.0105; T y3 : PO y 

T a : ISE = 1.0; T l2 : ISE = 1.2; T l3 : ISE = 1.4 



5.0 h 10 
0.020 



Figure 2. Tuning chart of systems without time delay 



3. Systems with time delay 

3.1. Stability region. The stability zones are defined in [5] and [S] as follows: 
(1) The parameter h must be included in the range from zero to h u given by 

(3.1) h u = - cos(z a ) + t p z a sin(z a ) 

where z a is the first positive solution of 



tan(z a ) = — 



(2) The parameter hi, for a given value of h, must satisfy the following inequal- 
ities 

(3.2) S r (zi) < and8 r {z2) > and hi > 

where 

8 r (z) = hi — zsin(z) — t p z 2 cos(z) 
and Z\ and z 2 are the first two positive roots of 
h + cos(z) — t p zsin(z) = 
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The stability region borderline T s is plotted in Fig. 3 for t p = 0.55 and in Fig. 4 



for t p = 2.5. 



It is also convenient to consider the stability phase margin (PM), whose expres- 
sion, deduced from (jl.5p . is 



tan(z& + PM) 



hi + z\ht p 



z b (h - hitp) 
where z b is the solution of 

h 2 + % =i + ty b 

Z b 

The curves relative to PM = 30°, PM = 45° and PM = 60° are also included in 
Fig. 3 and Fig. 4. 

3.2. Tuning. The performance indexes PO y , PO v and ISE are evaluated in a 
period equal to seven times the process time delay, which is considered enough to 
obtain sound results, in order to reduce the computing complexity to an acceptable 
level. Therefore, for given values of h and hi, the controlled variable y(t n ) and 
afterwards the controller output u(t n ) multiplied by K, named v(t n ), are calculated 
in 701 equally spaced points from t — to t = 7 according to (|B.3[) and (IB. 41) . 
The opposite of the minimum negative value of y(t n ) is assumed as PO y and the 
analogous value of v(t n ) as PO v . The ISE is computed with the triangle rule by 
means of (IB. 51) . In both Fig. 3, valid for t p = 0.55, and Fig. 4, valid for t p = 2.5, 
the curves T y for PO y — 0.0105 and T v for PO v = 0.1 are added. The optimum 
tuning point, named -B4, has been determined as follows: 

• let C y and C' v denote respectively each point lying on the curves T y and 
r„, corresponding to the 50 equally spaced values of h from h — to the 
maximum allowed by the stability 

• compute the ISE in C y , if T y is below or in C v , if above, since both 
PO y and PO v increase when, for a given h, hi increases 

• B4 is the point having the minimum value of the ISE 

It has been found that for t p — 0.55 (Fig. 3) B4 lies on the curve T y and for t p = 2.5 
(Fig. 4) on r„. 

3.3. Comparison. The proposed optimum tuning is compared with the Ziegler- 
Nichols time-domain and frequency domain [5] and with Zhuang-Atherton ISTE 
for setpoint change one [6]. Taking into consideration the equalities a = K/t p , 
K U K = (1 + (z u t p ) 2 ) ' 5 , T u /L — 2ir/z u , the parameters are given by: 

(1) Ziegler-Nichols time-domain 

0.9 _ 0.9 
~a ~ ~K 

(2) Ziegler-Nichols frequency-domain 



K p = — =^rt p ; Ti/L = 3 



K p = 0AK u = ^ (1 + zlt 2 p )°- 5 ; TJL = 0.8TJL = 0.8 2ir/z u 

where z u is the first positive root of tan(z u ) = —z u t p 
(3) Zhuang-Atherton ISTE setpoint 
(a) 0.5 <= tp <= 0.9 

K p = (1/tp)- - 559 ; Ti/L = t p /(0.883 - 0.158 ftp) 
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FIGURE 3. Tuning chart of systems with time delay- t p = 0.55 



(b) 1 <= tp <= 10 
71 2 

K p = — — (l/g-°- 921 ; Ti/L = i p / (0.968 - 0.247/t p ) 

A 

Since h = K K p and hi = hL/Ti, it is possible to calculate the optimum tuning 
points B\, B 2 and B 3 , also plotted in Fig. 3 and Fig. 4, for these three tuning 
rules. 

The optimum controllers parameters h and hi are listed in Table 1 whereas the 
performance indexes PO y , PO v and ISE are in Table 2 for some values of t p , from 
t p = 0.1 to t p = 10. From Table 2 it follows that: 

(1) the controlled variable response is overdamped if t p >= 0.85 for the pro- 
posed method and always for the others ones. 

(2) the controller output response is overdamped, if t p <= 1, for the two 
Ziegler- Nichols and Zhuang-Athcrton methods. Moreover the values of PO v 
for t p > 1 are unacceptable for the first two methods since higher than 0.10, 
but good for the third one since lower than 0.06. 

(3) neglecting the Zicgler-Nichols methods for t p > 1 in accordance with the 
previous consideration, the proposed method allows the best values of the 
ISE. 



A NEW OPTIMUM TUNING METHOD IN TIME-DELAY SYSTEMS 



7 




Figure 4. Tuning chart of systems with time delay - t p = 2.5 



Table 1. Optimum controllers parameters 



tp 


h 


hi 


Z-N 
time 


Z-N 
freq. 


Z-A 
ISTE 


Prop. 


Z-N 
time 


Z-N 
freq. 


Z-A 
ISTE 


Prop. 


0.10 


0.090 


0.416 




0.45 


0.030 


0.237 




0.787 


0.25 


0.225 


0.475 




0.50 


0.075 


0.243 




0.738 


0.40 


0.360 


0.552 




0.60 


0.120 


0.262 




0.724 


0.55 


0.495 


0.636 


0.563 


0.70 


0.165 


0.285 


0.609 


0.737 


0.70 


0.630 


0.724 


0.644 


0.92 


0.210 


0.311 


0.605 


0.763 


0.85 


0.765 


0.814 


0.718 


1.10 


0.255 


0.338 


0.589 


0.766 


1.00 


0.900 


0.905 


0.786 


1.15 


0.300 


0.365 


0.570 


0.744 


2.50 


2.250 


1.835 


1.656 


2.10 


0.750 


0.654 


0.576 


0.682 


4.00 


3.600 


2.774 


2.553 


3.00 


1.200 


0.947 


0.578 


0.654 


5.50 


4.950 


3.715 


3.423 


3.80 


1.650 


1.241 


0.574 


0.633 


7.00 


6.300 


4.656 


4.274 


4.75 


2.100 


1.535 


0.569 


0.628 


8.50 


7.650 


5.598 


5.111 


6.00 


2.550 


1.829 


0.564 


0.640 


10.00 


9.000 


6.540 


5.936 


6.65 


3.000 


2.123 


0.560 


0.622 



Using a least squares fit, the following empirical formulae, related to the proposed 
method, are obtained from the results shown in Table 1, namely 
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Table 2. Optimum controllers performance indexes 



tp 


POy 


PO v 


ISE 




Prop. 


Z-N 


Z-N 


Z-A 


Prop. 


Z-N 


Z-N 


Z-A 


Prop. 






time 


frcq. 


ISTE 




time 


freq. 


ISTE 




0.10 


0.010 








0.014 


6.095 


3.186 




1.524 


0.25 


0.010 








0.029 


5.192 


3.266 




1.674 


0.40 


0.010 








0.044 


4.605 


3.261 




1.788 


0.55 


0.010 








0.086 


4.193 


3.229 


1.998 


1.869 


0.70 


0.010 








0.099 


3.896 


3.191 


2.098 


1.945 


0.85 










0.100 


3.674 


3.154 


2.212 


2.037 


1.00 










0.100 


3.502 


3.119 


2.333 


2.129 


2.50 




0.177 


0.112 


0.032 


0.100 


2.822 


2.925 


3.089 


2.939 


4.00 




0.849 


0.609 


0.049 


0.100 


2.647 


2.863 


3.698 


3.582 


5.50 




1.523 


1.108 


0.053 


0.100 


2.574 


2.837 


4.181 


4.077 


7.00 




2.197 


1.608 


0.053 


0.100 


2.536 


2.825 


4.555 


4.458 


8.50 




2.872 


2.108 


0.052 


0.100 


2.513 


2.819 


4.849 


4.754 


10.00 




3.548 


2.608 


0.050 


0.100 


2.498 


2.815 


5.084 


4.993 



• optimum tuning points lying on T y for 0.1 <= t p <= 0.7 

h = 0.4541 - 0.1035 t p + 1.0794 t 2 p ; hi = 0.8271 - 0.4805 t p + 0.5613 1\ 

• optimum tuning points lying on T v for 0.85 <= t p <= 10 

h = 0.5884 + 0.5826 t p + 0.0033 t 2 p ; K = 0.7874 - 0.0434 t p + 0.0028 1\ 

The values of the parameters and performance indexes of the proposed method, as 
per these empirical formulae but considering all the decimal digits of the coefficients, 
are listed in Table 3. 

Table 3. Empirical parameters and indexes of the proposed method 



tp 


h 


hi 


PO y 


POy 


ISE 


0.10 


0.4546 


0.7846 


0.0096 


0.0100 


1.5259 


0.25 


0.4957 


0.7420 


0.0179 


0.0361 


1.6700 


0.40 


0.5854 


0.7247 


0.0103 


0.0512 


1.7864 


0.55 


0.7237 


0.7326 


0.0096 


0.0701 


1.8762 


0.70 


0.9106 


0.7657 


0.0104 


0.1052 


1.9410 


0.85 


1.0861 


0.7525 




0.0850 


2.0509 


1.00 


1.1744 


0.7468 




0.0974 


2.1315 


2.50 


2.0658 


0.6965 




0.1294 


2.8900 


4.00 


2.9722 


0.6589 




0.1115 


3.5606 


5.50 


3.8935 


0.6340 




0.0913 


4.0928 


7.00 


4.8298 


0.6218 




0.0821 


4.4885 


8.50 


5.7810 


0.6224 




0.0895 


4.7722 


10.00 


6.7473 


0.6357 




0.1158 


4.9695 
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4. Conclusion 

In this paper a new deadbeat response tuning method is applied to the PI con- 
trollers in first-order time-delay systems and compared with some existing ones. 
The involved performance indexes are strictly calculated by means of a new an- 
alytical solution of the differential difference equations relative to a step setpoint 
response. Moreover tuning charts, having as coordinates the two PI parameters 
and provided with the stability zone borderline and with the performance indexes 
curves, have been introduced. These charts, previously prepared, may be considered 
a general and an effective tool, which does not require any advanced or theoretical 
knowledge, for the start-up and the maintenance of any controller. 

Appendix A. Step setpoint response in systems without time delay 

During a step setpoint change from r = to r = 1 applied to the closed- 
loop system in steady condition, the differential equation, relative to the controlled 
variable y and deduced from ljl.4|L and its underdamped solution are given by 



(A.1) +U (l + h) m +h y(t) = h 

(A.2) y(t) = 1 + e- at (-cos(bt) - ~ sin(bt)) 

where 

a = 0.5(l + /i) 

(A.3) b=— (AhU~tj(l + h) 2 ) - 5 

H 

For the controller output u(t) multiplied by K, named v(t), it follows from (jl.lj) 
and (|A.2j) 



dt 

1 + e~ at (-cos{bt) + ~ a + L [ ' " sin(bt)) 



( A - 4 ) 



b 

The derivatives of y(t) and v(t) with respect to t are given respectively by 

(A.5) « = e-«^ «n(6t) 

at b 

(A.6) ^ =e- at (a 2 + b 2 )(cos(bt)-^^ sin(bt)) 

dt b 

Introducing the solution t of (|A.5[) . equated to zero, in (|A.2|) . we obtain for the 

overshoot PO y 

(A.7) PO y = e-™/ b 

Analogously introducing the solution t of (|A.6j) . equated to zero, in (|A.4j) . we obtain 
for the overshoot PO v 

(A.8) PO v = e~ {a ^ arctan(6/(a-l)) (1 _ 2 fl + O 2 + 6 2 ) ' 5 

From (|A.2| it follows 

0.25 5a 2 + 6 2 



(A.9) ISE= / {y{t)-l) 2 dt 



b 2 
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Appendix B. Step setpoint response in systems with time delay 

During a step setpoint change from r = 1 to r = applied to the closed- 
loop system in steady condition, the differential equation, relative to the controlled 
variable y and deduced from (| 1 . 5[) . is given by 

dy(t) +t d 2 y(t) 

= (hi+h^-)(r(t-l)-y(t-l)) fort>=0 
at 

subject to an initial condition of the form 

(B.2) y(t)=yo(t) = l for -l<=t<=0 

The analytical solution of (|B.1|) . obtained with the method of steps in [T], consists 
of a set of functions y n (t n ), each valid for n — 1 < t < n, whose independent variable 
t n has the temporal origin at t = n — 1, given by 

/-n o\ i=n— 1 j=n—2 

Vn{t n )= + e~ Wt > n •■■!<>■ f° r n>1 

In Example no. 2 of [T] there are also the recursive expressions suitable for the 
evaluation of the coefficients A n> i and B n j. 

For the controller output u n (t n ) multiplied by K, named v n (t n ), it follows from 
CCD) 

(B.4) v n (t n ) = y n+1 (t n+1 ) + tp d Vn+iitn + i) 

U.tn+1 

Finally the integral of the squared error ISE is calculated with 

ISE = + ~ o.oi ( yi (o)) 2 - \ o-oi Mi)) 2 

(B.5) ,i=7t=100 

+ °- 01 E E iVni^r)) 2 

r — 1 
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